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Abstract 
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case of an elliptic-cylindrical boundary. In the case of massless vector fields, an exact solution 
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energy of the bag should stabilize the bag against deformations from a spherical shape. 
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I. INTRODUCTION 



The Casimir effect is a manifestation of the quantum nature of fields. The presence of 
physical boundaries gives rise to changes in the zero point energy of a quantum field 
This has been studied theoretically for a number of geometries |2|, y], however, apart from 
fairly simple cases such as parallel plates, spheres and cylinders, exact results are difficult 
to find because analytic solutions to the field equations are required. The possibility of 
technological applications of quantum vacuum fluctuations |5| has made further inves- 
tigation of the Casimir effect in situations where the boundaries are either not perfectly 
conducting or of lower symmetry timely. 

The problem of the Casimir self-stress on a spherical shell was first proposed by Casimir 
611 as a possible mechanism for stabilizing a semiclassical model of the electron. However, 

n 

Boyer showed [7| that the zero point energy is positive in this case. While the concept 
of the zero point energy of a single object must be carefully defined to avoid ambiguity, 
this notion has been applied to Kaluza-Klein theories js, ol and cosmological models, for 



example in 
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llj . Another area where the idea of the zero-point energy of a single 



object has found application is in hadronic physics. There have been many papers on the 
zero-point energy of coloured fields confined in a spherical cavity following the early work 
of Bender and Hays 12| and Milton pj|]. Also, zero-point energy concepts have been 
applied to soliton 1J| , hybrid bag 15[ and string models pL6l] of hadrons. 
In recent years there has also been experimental and theoretical progress on the Casimir 
brce between separated objects. We refer the reader to review articles by Lamoreaux 
5|, Milton 171] and Klimchitskaya, Mohideen and Mostepanenko ISj. Of particular note 
is the recent experimental finding [loi] that the Casimir or, more properly, the Casimir - 
Lifshitz force can be repulsive, confirming a prediction of Lifshitz and collaborators {20 1 
for a system where the two plates have different dielectric permittivities and the vacuum 
is replaced by a liquid with a permittivity whose value is between those of the plates. 
It is well known that the Casimir energy is sensitive to geometry - even the sign of the 
energy is difficult to predict. It is natural, therefore, to investigate the modification of zero 
point energy when the boundary is modified. Experimental investigations 21[ have shown 
that changing the boundary on the nanometre scale alters the Casimir - Lifshitz force by 
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a factor different from that predicted usin^ tfie proximity force approximation, though 
more sophisticated calculations 



22 



work 



23| can account for the change in the force. Recent 



24| has provided experimental confirmation of an exact theoretical calculation of 



the effect of varying the amplitude of a sinusoidally corrugated surface. 
Kitson and Romeo 28|] looked at the theoretical problem of an infinite cylinder de- 
formed from circular to elliptical cross-section. For small ellipticity (or eccentricity) e 
a perturbative-like expansion for the zero point energy can be found: 
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Sc{a,e) = Sc{a,0) 



1 + -e^ + 0{e') 
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where a is the semi-major axis of the ellipse, and Sc{(i,0) is the regularized zero point 
energy per unit length of an infinite circular cylinder of radius a, which has the numerical 



value 



29| 



£ciO) 



0.01356 



(2) 



This agrees with earlier work of Kvitsinsky 



30l |. Noting that an ellipse with semi- major 



axis a and ellipticity e has area vra^-y/l — e^, Kitson and Romeo's result can be written in 
terms of R, the geometric mean of the semi-major and semi- minor axes R = \/ab 

0.01356 



S(e) 



[1 + 0{e 



(3) 



that is, if the cross-section of the cylinder remains constant under the deformation of the 
cross section from a circle of radius R to an ellipse, then the zero point energy per unit 
length does not change, up to C(e^). This is a remarkable result, and begs the question 
whether it can be generalized to other deformations of highly symmetric boundaries. 
In this paper we will consider the deformation of a sphere to a spheroid. It is known that 
the zero point energy for a spherical boundary is positive [7| , whereas it is negative in the 
cases of parallel plates and infinite circular and elliptical cylinder geometries. As these 
geometries are limiting cases of spheroidal geometry, it is interesting to investigate the 
change in zero point energy as a sphere is deformed into a spheroid, and to see whether 
this change lowers or increases the zero point energy. Some care may be needed here, 
as while the parallel plate configuration is the limiting case of an oblate spheroid with 
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ellipticity taken to unity, there is also a change in topology, so the value of a physical 
quantity such as the zero point energy may not be continuous in this limit. 
It should be recognized that the situation we are discussing is an idealization where 
boundary conditions have replaced a physical interaction, and care needs to be taken, 
especially with the appearance of divergences. An alternative approach is to replace 
the boundary with a renormalizable coupling between the fluctuating field and a non- 



dynamical smooth background field [25l-l27l| . This approach allows the use of the standard 
renormalization procedure of quantum field theory without boundaries to obtain finite 
results. 

In section 2 we consider the case of a massless scalar field subject to Dirichelet and Neu- 
mann boundary conditions on the surface of a prolate spheroid. We obtain the zero point 
energy as a perturbative expansion in the ellipticity. In section 3 we study the Maxwell 
equations with perfect conductor boundary conditions on the surface of a spheroid. Un- 
brtunately an analytic solution for the fields is not possible in this case. In previous work 
3l|, we proposed using a zonal approximation by weighting components of the known 
axial symmetric field solutions with a spheroidal boundary. However, we will show that 
this approximation is poor, and likely to give the wrong sign for the first correction term 
in the perturbative expansion of the zero point energy. In section 4 we examine the im- 
plications of our results for the MIT bag and flux tube models of hadron structure in 
QCD. 

II. MASSLESS SCALAR FIELD IN SPHEROIDAL GEOMETRY 

The spheroidal boundary is given by the Cartesian equation 

5+^-1 (4) 



a 



where both a and h are non-zero. A prolate spheroid has a> h whereas an oblate spheroid 
has a < b. We shall consider a prolate spheroid; the results for an oblate spheroid can 
be easily obtained from those we derive here. We can think of the prolate spheroidal 
boundary as the solid figure formed by rotating the ellipse with semi-major axis a and 
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semi-minor axis b about the The elhpticity of the spheroid is 



^ (5) 



a" 



It is most convenient to work in prolate-spheroidal coordinates 32|, which are 

related to Cartesian coordinates by 

X = fiv 

y = /V(e-l)(l-^')cos0 

z = /v/(e-l)(l-r/2)sin0 (6) 

where / = \/a? — IP' = ea. The domains of the prolate spheroidal coordinates are 1 < 
^<oo,— 1<?7<1 and < cj) < 2n. In these coordinates the boundary @ is given by 

5 = i. (7) 

Let us consider a real massless scalar field (f that satisfies the homogeneous Dirichelet 
condition on the boundary of the prolate spheroid. In natural units the field equation 
and boundary condition are 

{V^ + u;^)ip = (8) 

^k=l/e = (9) 

where u is the frequency. In prolate spheroidal coordinates the field equation is separable, 
and the general expression for the field in the interior of the boundary can be found 33 1 



oo I 
/=0 m=-l 

where 7 = fu. Here the functions S^^^^ are the radial prolate spheroidal functions of the 
first kind and ps™^^"* are the angular prolate spheroidal functions of the first kind 32, 34 1. 
The homogeneous Dirichelet boundary condition is satisfied when 

^™{l)(i,;,V)=0 (11) 

where we have introduced the dimensionless frequency z = ua. Exterior to the boundary 

the general expression for the field is similar: 

00 I 

V'"'°(e,r],0) = E E Brsr^'\^,j')psf'\v,j')exp{zm<P) (12) 



where Si are the radial prolate spheroidal functions of the third kind. 
With these solutions to the field equation we can proceed to finding the zero point energy. 
We will use two methods to find the zero point energy: the zeta function method and the 
Green's function method. 



Zeta Function Method 



The zeta function for the field is given by 



(13) 



where are the eigenenergies of the field, s is a complex variable and /i is a positive 
mass scale. The zero point energy is related to the principal part of the zeta function 



^C=^PPs=-l/iC(s) 



(14) 



Let the positive solutions of the spheroidal boundary condition equation f lTT]) be given by 
zimn- For large enough Re(s) we can write the zeta function for the interior field as 

oo I oo 



C''%,a,e) = {f,arYl E E 



Imn 



(15) 



/=0 m=—l n=l 

The summation over n can be performed via the argument principle 
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n=l ^ 



Sr'^'^'il/e, zV) 



s: 



m(l). 



(16) 



Here the prime denotes differentiation with respect to z, and the contour is a closed path 
that encloses all the positive zeroes of the radial spheroidal function and avoids the origin. 
We can make progress via a formal asymptotic expansion of the radial prolate spheroidal 
functions for small ellipticity in terms of spherical Bessel functions of the first kind ji{z) 

+ l + m^ -1 



Sr^'\l/e,z'e') ~ Uz)- 



4/2 + 4/ - 3 

't^ + 1 — 3m^ m 

8P + 81-6 + y ' ^'^^^ 



(17) 



6 



The integrand of the zeta function can now be expressed using this expansion as 



ji{z) 4/2 + 4/ -3 
The leading order term of flTSi 
spherical boundary of radius a 



SZ '-i^ + z 



Ji[z) V Ji[z) 
can be recognized as the integrand we would obtain for a 
. The O(e^) term has two terms: the first is proportional 



to the leading order term, the second is continuous on C and has an antiderivative, so 
will integrate to zero. We can perform the summation over m as the leading order term 
is independent of m, and for the next-to-leading order term we can use the relation 

>^ + l + _ 2l + l 

^ 4/2 + 4/-3 ~ 3 ■ ^ ^ 

m,=—l 

So now we have that the next-to-leading order term is simply —■^se^ times the leading 
order term, and we can write the zeta function as 

C^'°(., a, e) ~ C''^(., a, 0) (l - ^e^ + 0(e^)) . (20) 

For the exterior field, we can do a similar analysis. The formal asymptotic expansion of 
the radial prolate spheroidal functions S^^^^ is the same as that for S^^^^ except that the 
spherical Bessel functions of the first kind are replaced by spherical Hankel functions of 
the first kind h\^\ Thus we obtain the expansion of the zeta function in the exterior of 
the spheroidal boundary: 

C"'°(s, a, e) ~ C"'°(s, a, 0) (l - ^e^ + 0{e')) . (21) 

The total zeta function is the sum of internal and external contributions 

C°(s,a,e) = C''°(s,a,e) + C"'°(s,a,e) 

~ C°(s,a,0)(l- V + 0(e^)). (22) 

The first term is the total zeta function for a massless scalar field satisfying the ho- 
mogeneous Dirichelet condition on the boundary of a sphere of radius a. This term is 
well-behaved at s = — 1, and taking the principal part gives the zero-point energy 



£^(a, e) ~ f ^(a, 0) ( 1 + -e^ + O(e^) ) (23) 



where the leading term is the zero point energy of the massless scalar field for a sphere 
with Dirichelet boundary conditions 37 1 



^n. ^ 0.00281... , 

S^{a,0) = (24) 

a 

We can also consider Neumann (or Robin) boundary conditions. The homogeneous Neu- 
mann boundary condition for the massless scalar field in prolate spheroidal geometry will 
be satisfied when the radial function satisfies 

^Sr^'\^,u'e%=^/^ = Sr{u,e) = 0, (25) 

where u is now the dimensionless frequency. We can find the zeta function for the interior 
region 

oo I oo 

C^-^(.,e) = (H^^EE<- (26) 

1=0 m=—l n=l 

As before, we can convert the summation over n to a contour integral over complex u 

Here the prime denotes differentiation with respect to u. Again using the formal expan- 
sion of the radial prolate spheroidal functions for small e, equation (fT7|) . we obtain an 
expression for the integrand of the zeta function which is very similar to equation f|T8|) . 
except that the Bessel function ji is replaced by its derivative, the derivative j'l is replaced 
by the second derivative, and we find an extra term which is a power of u: 

u — ~ 

Sr(n,e) 



^ S 1 ^ SJI \ / 



j'liu) 4/2 + 4/ - 3 



su ,, , + u 



As the contour avoids the origin, the u"''"^ term will integrate to zero, and the analysis 
can proceed as before: the leading term being the same as for a spherical boundary, and 
the next-to leading-term having one part which integrates to zero and another part which 
is proportional to — |s times the leading term. Hence we arrive at 

C^-^(., a, e) ~ C^'^(., a, 0) (l - ^e^ + 0{e')) . (29) 



Clearly this analysis extends to the field exterior to the prolate spheroid, and so we see 
that equation fl23|) holds for both Dirichelet and Neumann boundary conditions and can 
be considered a general result for the zero point energy of the massless scalar field. 



B. Green's function method 



The zero point energy is related to the reduced Green's function in a region of space 

= J— J dco j d^xuj^g{x,x). (30) 
For the interior of the spheroidal boundary, with the homogeneous Dirichelet boundary 



condition, we have the reduced Green's function 



g^'^(x,x') = 

oo I 
1=0 m=-l 



where 



2/ + l(/-m)! m(i) 2n /• 

-pS; (?7, 7 j exp[tm(, 



(32) 



Att {l + m)V 

For the exterior region, the reduced Green's function is similar, except that the radial 
prolate spheroidal functions of the first and third kind are interchanged. The spatial 
integral in equation (l30ll becomes 



I(..e)^/d3..r(x.x 

OO I 

--EE 



1=0 m=-l 



dvXrivA,z'e')Xr{vA,z'e') 



-1 



X 



Sf^\l/e,z'e^) 



dae-v'z'e'){sr^'\^,z'e') 



(33) 



S;'^'\l/e,z^e^) 



ze 



ze 



where we have rescaled ^ by a factor of ze. As in the previous subsection, we now make 
formal expansion of the radial and angular prolate spheroidal functions. Note that the 
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argument of the spherical Bessel functions in the expansion fll7l) becomes ^. For the 
angular prolate spheroidal functions we have the expansion for small 7^ 

We obtain the expansion of the spatial integral for small ellipticity 

X(2,e) ~X(°)(-2)+2:^'H^)e' + C(e^). (35) 
As we expect, the leading term is the same as the spatial integral for the same problem 



on the boundary of a sphere of radius a [39| : 

00 I „2n 



X 



' duu {ji{u)) + / duu \^hl'{u)j 



(36) 



Ji{z) Jo h\'>{z) 
The next-to leading order term is given by 

X(^)(z) = limi|^X(^,e). (37) 

We note that in the integral over the angular variables, only the leading term contributes 
because of the orthogonality of the associated Legendre functions P[". For the radial 
integrals, we need to use Leibnitz's rule. There is one term involving r/^, however this can 
be dealt with using the recurrence relation for associated Legendre functions 

(2/ + l)vPriv) = {l + l- m)PllM + (/ + m)Pl1M. (38) 
and the orthogonality relation. After further simplification we can write 

X(2)(.) = lx(°)(.) - 1^ (zX^'\z) + QT{z)\ . (39) 

When the second term is integrated over z, we can use an ie prescription to deform the 
contour in the complex plane to avoid the poles along the real axis, and close the contour 
around the upper half-plane. The integral will them vanish as it is a total derivative. 
Thus the next-to-leading order term is times the leading order term, which agrees 
with the above calculation using the zeta function method. 
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C. Discussion 



Our equation ( 123|) for the zero point energy of the prolate spheroid is a new result. We 
can generalize this result to an oblate spheroid. Rather than repeat our derivation above, 
we note that the ellipticity e' of an oblate spheroid is related to that of a prolate spheroid 
with a ^ b hj 

Substituting into equation (l23l) yields 

E°(a, e') ~ E°(a, 0) (^1 - + 0(e^)) , (41) 
similar to the result for a prolate spheroid. 

We can gain more insight into our result by considering a sphere of radius R. The volume 
of this sphere is equal to that of the prolate spheroid when 

R = a^l-e2, (42) 

which gives the zero point energy of the spheroid as 

£(a,e) ~£(i?,0) (l + 0(e^)), (43) 

i. e. the leading correction term vanishes when the sphere is deformed to a spheroid in 
a manner that preserves the volume. This agrees with the result of Kitson and Romeo 



for the elliptical cylinder geometry 



28 



It is natural to speculate that this result may 
hold more generally |41j], though this is difficult to investigate. In the case of the elliptical 
cylinder, the result could be derived using a conformal map, which only requires continuity 
of the field solutions and is independent of boundary condition. For the spheroid we do 
not know the range of validity for our expansion in the ellipticity, so confirmation of this 
result may require experimental input. 

In the limit e — )■ 1, a prolate spheroid becomes a circular cylinder, which has negative zero 
point energy. We note that our result for the zero point energy of the prolate spheroid is 
positive and increases in magnitude as the ellipticity increases. For an oblate spheroid, 
which becomes parallel plates in the limit e' — )■ 1, our result for the zero point energy does 
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decrease as the ellipticity increases, however, it remains positive. (A plot in an earher 



paper 



3l| erroneously showed a negative zero point energy.) So the change in sign of the 
zero point as the geometry changes must come from higher order terms in the ellipticity 
expansion, or be associated with the change in topology at the limit. In figure 1 we show 
the total zero-point energy for a scalar field subject to Dirichelet boundary conditions on 
both a prolate and an oblate spheroid. 

III. ELECTROMAGNETIC FIELD IN SPHEROIDAL GEOMETRY 



In general Maxwell's equations are not separable in spheroidal coordinates 33|. This 
means that we cannot use either zeta function or Green's function methods to investi- 
gate the zero point energy of the electromagnetic field subject to appropriate boundary 
conditions on the boundary of either a prolate or oblate spheroid. However, there is one 
known case of the Maxwell equations separating in prolate spheroidal coordinates, which 



is when the field is axially symmetric 40|]. In a previous paper we proposed that the 
zero point energy could be obtained in a zonal approximation by suitably weighting the 
zeta function of the axially symmetric solution, and obtained a result with the leading 
correction term in the ellipticity expansion having opposite sign to the spherical result 



3l| . We now show that this result is unreliable, and that the zonal approximation is 
probably a poor conjecture. 

In the case of a scalar field we can compare the exact result with the result of the zonal 
approximation. The zeta function for the field satisfying the homogeneous Dirichelet 
boundary condition on the prolate spheroid is given by 

The zonal approximation is made by setting m = in the integrand and replacing the 
sum over m by a factor (2/ + 1). We denote the zonal approximation to the zeta function 
by (, and use the expansion in ellipticity, equation (fT7|) to write 

({s, a, e) ~ C^°)(., a) + C^'\s, a)e^ + 0{e'). (45) 

For simplicity we set the mass scale /i equal to a, and lei v = 1 + 1/2. The leading order 
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term is the zeta function for a spherical boundary 



c 



(0), 



2z/ 



s, a 



1=0 



dz z 



c 



(46) 



which has a Laurent expansion about s = —1 42 1 

1 1 



3157rs + 1 

The next-to-leading order term is found to be 

5 - 4z/2 2z/ 



0.00889... + C(s + 1) 



(47) 



s, a = — s 



E 



16(1 - z/2) 27ri 



c 



(45 



This expression needs analytic continuation from the domain where it is absolutely con- 
vergent Re(s) > 2 to — 1 < Re(s) < 0. This can be done by considering the contour to run 
along the imaginary axis from z = +iy = +ioo to z = —iy = —ioo, avoiding the origin 
via a semi-circle of radius e, and being completed by a semi-circle in the right half-plane, 
see figure 2. The line integrals along the large semi-circle, radius p oo and the small 
semi-circle, radius e both vanish. The integral along the imaginary axis becomes 



c 



(2)/ 



^8(1 



1=0 



^_sm — / dyyU-^^-l + — 



(49) 



where ii{z) are the modified spherical Bessel functions. Now continuing to a neighborhood 
of s = — 1 we find the Laurent expansion 

3 1 2561 - 18907- 5670 In 2 1 



- 0.03421 . . . + C(s + 1) (50) 



647r (s + 1)2 403207r s + 1 

where 7 = 0.57721 ... is the Euler-Mascheroni constant. Taking the principal part of 
both the leading and next-to-leading terms, and factoring out the result for the spherical 
boundary, we have the zero point energy using the zonal approximation 



^^•°(a, e) ~ ^^'°(a, 0) (l - 3.85312 . . . + C(e^)) . 

However, we can expand the zeta function, equation (jS]), for small ellipticity as 

00 I 



(51) 



00 t / 

C^-°(.,a,e)~5^ J]C(^,«,0)(l 

/=0 m=-l ^ 



4P + 41-3 



(52) 
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where C;^^ (5,0,0) is the term appearing under the summation in equation ( H^ for the 
zeta function for the spherical boundary. The summation over m is performed using 
equation (fT9l) giving 

oo 

C^-°(., a, e) ~ 5^(2/ + l)C{s, a, 0) (l - '-e' + 0{e')) (53) 

1=0 

as expected from our earher result, equation fl20|) . Using the zeta function for the sphere, 
equation (HTIl we get the exact result for the zero point energy 

£^'^(a, e) ~ £^'°(a, 0) (l + 0.25759 ...e^ + 0{e^)) . (54) 

We see that the zonal approximation gives a next-to-leading term that is the wrong sign 
and an order of magnitude too large. Additionally, the Laurent expansion in the zonal 
approximation, equation (!50|) . has a double pole, which contradicts the Poisson kernel 
method. We thus conclude that the zonal approximation is a poor approximation in this 
case. We also note that if we had used the boundary deformation conjecture, equation 
( 123|) . which takes into account the external field, the next-to-leading order term would be 
+^S^'^{a, 0)e^ which has the correct sign, and is not too different in magnitude from the 
exact result. 

In the coordinate systems where Maxwell's equations are separable (e. g. spherical polar 
or elliptical cylindrical coordinates), the field equations can be broken down into two sets 
of equations with solutions corresponding to TE or TM modes. For perfectly conducting 
boundary conditions, each of these modes (up to polarization states) is equivalent to 
a massless scalar field subject to either homogeneous Dirichelet or Neumann boundary 
conditions. So, in these geometries, finding the zero point energy of the massless scalar 
field with either Dirichelet or Neumann boundary conditions gives the zero point energy 
of the electromagnetic field as the sum of the zero point energies for the two solutions of 
the massless scalar field. In the case of the spheroid, we were able to find the zero point 
energy for both solutions of the massless scalar field, so it is interesting to conjecture 
that, for small ellipticity, the zero point energy of the electromagnetic field will follow our 
boundary deformation conjecture 

^™(a, e) ~ f ™(a, 0) (^1 + ^e' + 0(e^)^ (55) 
14 



where S^^{a,0) is Boyer's result for a sphere of radius a 

^,^^ 0.04617... 

a 

which is the sum of zero point energies for Dirichelet and Neumann boundary conditions 
minus the contributions from n = modes. The difficulty with proving this conjecture 
is that we do not know exactly what form the boundary conditions will take for the 
solutions of the field equations in either prolate or oblate spheroidal coordinates, as these 
solutions do not separate into 'radial' and 'angular' parts. However, the simplicity of 
the perfectly conducting boundary conditions, which require only the 'radial' part of the 
solution to vanish, or have vanishing first derivative, on the surface gives us some grounds 
for believing that these boundary conditions will not differ greatly from those we have 
used for the massless scalar field, and so the boundary deformation conjecture should be 
applicable for the electromagnetic field. 

IV. ZERO POINT ENERGY IN THE MIT BAG MODEL 



The MIT bag model 43| is a phenomenologically successful model of hadron structure. 
In the model, the QCD vacuum has infinite colour magnetic permeability fi, while inside 
the bag the permeability is unity. This confines colour magnetic and electric fields to the 
interior with boundary conditions on the bag surface 

n-E = = nxB (57) 

where n is a unit vector pointing in the inward normal direction. The original fit to 



hadron masses [4^ included a term —Z/R (where r is the radius of a spherical bag) to 
account for zero point energy and centre of mass effects, with a value of Z ^ 1.8 giving 
a good fit. We note that the fit for the strong coupling constant inside the bag typically 
gives fairly large values as ~ 2.2 in order that the colour hyperfine splitting reproduce 
the nucleon-delta mass splitting of 300 MeV. 

The zero-point energy of the confined colour fields, with divergences subtracted, can be 
estimated [3] 

S{R) ^ +^ (58) 
15 



which is of opposite sign to the value found in hadron mass fits. As the leading divergence 
is associated with the bag surface, more realistic boundary conditions incorporating a 



'skin depth', will not greatly alter this result 



The contribution of zero point energy 



of fermion (quark and antiquark) fields is small, and is usually neglected. Centre of mass 



motion is also estimated to give a positive contribution to hadron masses 45|, |46l]. One 
possible resolution to this is to include further terms in the mass formula, such as a 
constant force term FR ^ , which allows reasonable fits to the hadron masses. 
Another approach, in light of the large values of as found in the mass fits, is to consider 



non-perturbative effects. Recently Oxman, Svaiter and Amaral |48| considered modifying 
the gluonic Green's function to include confinement effects. Based on Schwinger-Dyson 
analysis, they introduced a modified reduced Green's function 

where A is a low momentum scale A ~ Aqcd and A is a positive constant, typically A ~ 
1. This modified Green's function reproduces the expected large momentum behaviour 
gik"^) — —l/k"^, and has a power law behaviour at low momentum ^(/c^) ~ {—k"^)^. The 
zero point energy can be expressed in terms of the reduced Green's function J2|] 

^ = -\Y.hl duj\ng{uj'-ul). (60) 

Now expressing the modified reduced Green's function in terms of the unmodified reduced 
Green's function 

g{e)=g{k^)-^g{k^-AY^^ (61) 

we can find the modified zero point energy 

8 = -A^(0) + (A + l)S{k^) (62) 

where £{m?) is the zero point energy calculated using the unmodified Green's function for 
a field of mass m. The second term will be small, so taking A ~ 1 we obtain the modified 
result 

^ ~ (63) 
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which is in better agreement with the usual hadron mass fits. 

Using the boundary deformation conjecture, we can investigate the stabihty of the spher- 
ical bag against deformations. We note that the zero point energy of the bag no longer 
includes terms coming from the exterior of the boundary, so the boundary deformation 
conjecture is unlikely to be exact. However, as we saw in the case of the massless scalar 
field in the interior region, the conjecture is able to give the correct sign and magnitude 
of next-to-leading terms in the ellipticity expansion. Let us first consider a meson state 
(quark - antiquark pair). At small separations of the quark and antiquark, lattice QCD 



simulations 49|-l5l| show that the pair sit in a roughly spherical potential, whereas at large 
separations the sphere becomes stretched into a fiux tube of roughly the same radius as 
the original 'bag'. We can model this as a spherical bag being deformed to a prolate 
spheroid with the semi-minor axis b staying approximately constant. Using the modified 
Green's function result above, we have the zero point energy of the spheroid 



S^-—ll+e' + 0{e^) . (64) 
a \ 6 J 

Now the semi-minor and semi-major axes of the spheroid are related by \? = a^(l — e^), 
so the zero point energy becomes 

0.7 1 



E{e) ~ (^1 - + 0{e')^ , (65) 

so increasing the ellipticity will increase the energy of the bag. We can also consider a 
baryon (three quark) state. At smal 



inter-quark separations, the quarks sit in a spherical 
bag. Lattice QCD simulations [5^ show that when one quark is separated from the 
other two a fiux tube is formed, similar to the quark -antiquark case, with the fiux tube 
again having approximately the same radius as the original bag. We can model this as a 
spherical bag being deformed to an oblate spheroid with constant semi-minor axis a. The 
zero point energy of the oblate spheroidal bag using the modified Green's function result 
is 

^ - (i - Y + ^^""0 ^66) 

where e' is the ellipticity of the oblate spheroid. Again we see that increasing the ellipticity 
increases the energy of the bag. These results are consistent with the fiux tube model 
result that the fiux tube has constant energy per unit length. 
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V. SUMMARY 



In this paper we have investigated the zero point energy of a spheroid. For a massless 
scalar field, subject to either Dirichelet or Neumann boundary conditions on the boundary 
of the spheroid, we were able to make an expansion in the ellipticity around e = and 
find the lowest order correction term for deviations from a perfect sphere. We used both 
zeta function and Green's function methods to obtain the lowest order correction, and 
both methods agree. Our result is consistent with the boundary deformation conjecture 



of Kitson and Romeo 



281], first seen in the case of an elliptical cylinder. In this case, the 



boundary deformation conjecture says that the zero point energy of a sphere does not 
change, up to 0{e^), when the sphere is deformed to a spheroid. In the case of a massless 
vector field (or electromagnetism), we were able to show that the zonal approximation 
we had previously used to calculate the zero point energy [3l| is likely to be incorrect, 
and we argued that there are good reasons for believing that the boundary deformation 
conjecture should apply in this case, even if an exact solution is not possible. Finally, we 
applied the boundary deformation conjecture to the MIT bag model, where we argued 
that zero point energy should increase when the bag is deformed from a sphere. 
It is tempting to relate the boundary deformation conjecture to the results of Ambj0rn 
and Wolfram 53|] for a cubiodal box boundary. However, the results for the cuboid are 
only calculated for the field on the interior of the boundary, and are thus both divergent 
and dependent on the renomalization method used (zeta function regular izat ion in this 
case), whereas our result for the spheroid (and the result of Kitson and Romeo for the 
elliptic cylinder geometry) includes both internal and external fields and is finite (as the 
divergences cancel) and unambiguous. The finite part of Ambj0rn and Wolfram's result 
shows that the sign of the zero point energy can change as the box is deformed, so it 
would be of considerable interest if their calculation could be extended to the external 
field, which would allow an unambiguous result. 

It would also be of interest to extend the work of this paper to spinor fields and fields 
with mass. Finally, the 0(e^) corrections in the ellipticity expansion should be calculated 
to enable us to place better limits on the boundary deformation conjecture. 
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Appendix: Spheroidal Functions 



The scalar Helmholtz equation is separable in prolate spheroidal and oblate spheroidal co- 
ordinates. As oblate spheroidal functions are easily related to prolate spheroidal functions 



32| . we shall only consider prolate spheroidal functions. The solution of the Helmholtz 



equation in prolate spheroidal coordinates has the form 



where the functions S and H satisfy the ordinary differential equations 



(67) 



d 



m 



d^ 



■d^ 
d 



.2 " 



m 



V 



Hiv) 




0. 



(68) 
(69) 



Here 7 = fu, m is an integer and \im is a constant of separation. We note that these 
equations are identical, however, they apply over different ranges of the variables, ^ > 1 
and |?7| < 1. The solutions to (!68l) are the radial prolate spheroidal functions; the solutions 
to ( l69l) that are regular and well behaved at 77 = ±1 are the angular spheroidal functions 
of the first kind 1 1"^) 



Angular Prolate Spheroidal Functions 

Angular prolate spheroidal functions of the first kind can be expressed as a series of 
associated Legendre polynomials 

00 

v^f'\v.i') = E i-^f<k{l')PJX2M (70) 

fe=— 00 
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where the series coefficients satisfy a three term recurrence relation 3J, |54| 



(71) 



with coefficients 



.m. 2^ ^ , (/-m + 2fc-l)(/-m + 2fc) 
''^^^ ^ ^ (2m + 4A;-3)(2m + 4A;- 1) 

R™r^2^ r/^9^V/^9^^n o.,2 (^ + 2A:)(/ + 2fc + l) + m^-l 

5,,(7 ) = (/ + 2fc)(Z + 2A; + 1) - 27 (2m + 4A:-l)(2m + 4A:-3) 

.m/ 2n 2(^ + "^ + 2/i; + l)(/ + m + 2A; + 2) 



-7 



(72) 



(2m + 4A; + 3) (2m + 4A; + 5) 

There are a number of normahzations for the angular prolate spheroidal functions: we 
have used the normalization of reference 34 1 



^ A m{l)/ 2\ m{l)( 2^ ^ (/ + m)! 

dr/ps; ' '(r/,7 )ps;/ ^(7^,7) = — -5w. 



(73) 



2/ + l(/-m)! 

In the spherical limit (7^ — > 0), equation (l69i) becomes the associate Legendre equation 
with eigenvalue A^m = ^(^ + !)• This enables us to make an expansion for the angular 
prolate spheroidal functions and their eigenvalues in powers of 7^ 



psr^^)(r^,7^) = Pr{ri)+l'Qfl + 0{i') 
W) = /(/ + l)+fA!S + (^(7')- 



(74) 
(75) 



Substituting these expansions back into the angular prolate spheroidal differential equa- 
tion fl69|) . and using the expansion in terms of associated Legendre polynomials (1701) it is 
straightforward, though tedious, to find 



mil)/ 2 



{l + m-l){l + m) 
2(2/- 1)2(2/ + 1) 
{l~m + l){l-m + 2) ■ 



f + 0(7^). (76) 



2(2/ + l)(2/ + 3)^ 

In figure 3 we plot the ratio of our expansion up to 0(7^) over the full angular prolate 
spheroidal function ps°''^''(r7 = 0.5,7^) for / = 0, 10,20 and 30 over the range < 7^ < 1. 
We see that the expansion is accurate to better than 1% over this range of 7. 



20 



Radial Prolate Spheroidal Functions 



The solutions to the radial prolate spheroidal equation ( 168|) are related to spherical Bessel 
functions. If we take the limit 7^ — )■ 0, ^ — t- 00 with 2; = 7^ = constant, and also make 
the substitution E{^) = (1 - we find 



z 



dz'^ dz 



+ 2z-f + (z2_A)/(^) = 



(77) 



which is satisfied by the spherical Bessel functions ji{z) and yi{z). We are interested in the 
radial spheroidal functions of the first and third kinds, which satisfy the physical boundary 
conditions at the surface of the spheroid and at infinity. These can be expressed as series 
of either spherical Bessel functions of the first kind ji{z) or spherical Bessel functions of 
the third kind hi{z) = ji{z) + iyi{z), also known as spherical Hankel functions of the first 
kind: [M] 



!■ y I k=—oo 

where f = j, h respectively, and 

00 

k=—oo 



(7J 



(79) 



The coefficients a]^ are the same as for the expansion of the angular prolate spheroidal 
functions. The normalization has been chosen such that for large 7^ 



Now we are able to make an asymptotic expansion in the ellipticity e = 1/^ 



^0) 



Sr^'\-^,z'e')r^Mz) + eX^:^ + 0{e^ 



.2^(2) 



(81) 



and then by substituting back into the radial spheroidal equation (168!) and using the 
Bessel function expansion fl78|) we can obtain 



Sr\l/e,z'e') ~ Mz)- 



't^ + 1 — 3m^ m 



+ 8/ - 6 2 



e^ + 0{e^). (82) 
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For the radial prolate spheroidal functions of the third kind Si (1/e, z^e^) we obtain the 
same expansion except that the spherical Bessel functions ji{z) are replaced by spherical 
Hankcl functions of the first kind. In figure 4 we plot the ratio of our expansion up 
to 0{e^) over the full radial prolate spheroidal function S'^'^^\\/ e, z^e^) for 2; = 10 and 
/ = 0, 10, 20 and 30 over the range < e < 0.3. We see that the expansion is accurate to 
better than 1% for e < 0.1, but can become much less accurate as the ellipticity increases. 
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FIG. 1: Total zero-point energy of a scalar field with Dirichelet conditions on a prolate spheroid 
(left) and an oblate spheroid (right) of small ellipticity. 
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FIG. 2: The contour of integration for the zeta function in equation 
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FIG. 3: The ratio of the approximation, equation (I76p . to the full angular prolate spheroidal 
function psl^^\ri = 0.5,7^) for Z = (full curve), 10 (dashed curve), 20 (dotted curve) and 30 
(dot-dash curve). 
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FIG. 4: The ratio of the approximation, equation (I82p . to the full radial prolate spheroidal 
function S^^^\l/e,z'^e'^) for z = 10 and / = (full curve), 10 (dashed curve), 20 (dotted curve) 
and 30 (dot-dash curve). 



28 



